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I .  INTRODUCTION 


A.  Background 

In  an  armory  or  a  depot,  in  transport  on  land  or  sea,  and  while  carried 
in  combat  vehicles,  munitions  are  packed  into  magazines,  pallets,  or  milvans. 
Such  logistic  stores  must  be  protected  against  the  risk  of  mass  detonation. 
Various  modes  of  interround  communication  of  violent  reaction  have  been  in¬ 
vestigated  by  Howe  and  co-workers1  b.  One  of  their  objectives  was  to  find  pre¬ 
dictive  criteria  which  could  be  used  to  reduce  the  vulnerability  of  munitions 
in  hazardous  environments. 

Generally,  munitions  are  sensitive  to  stimuli1 2  such  as  explosive  shock, 
air  blast,  fragment  impact,  and  fire.  An  accidential  explosion  of  a  single 
round  could  spread  throughout  the  entire  store.  This  is  referred  to  as  mu¬ 
nitions  fratricide  and  mass  detonation1’3 * * *.  Since  it  would  be  too  difficult 
to  assess  all  processes  involved  in  such  an  event,  we  consider  in  this  work 
the  shock- induced  phenomena  only.  V.e  have  chosen  to  do  so  for  two  additional 
reasons.  First,  shock  initiation  processes  dominate  the  early  stages  of  an 
interround  propagation  problem  especially  when  the  separation  is  small. 

Second,  these  processes  are  amenable  to  numerical  analysis  which  incorpo¬ 
rates  material  constants  and  data  in  reactive  hydrodynamic  codes.  We  have 
used  the  2DE  code7’8  to  carry  out  this  investigation  via  a  series  of  com¬ 
pilations.  As  experimental  data  are  scarce,  this  investigation  is  of  an 
exploratory  nature.  It  should  be  noted  that  our  calculated  results  furnish 
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considerable  details  which  are  not  readily  obtainable  using  other  approaches. 
These  results  also  disclose  certain  trends  which  may  be  used  for  predictive 
purposes . 

In  what  follows,  we  shall  formulate  the  munitions  sensitivity  problem  with 
input  data  adequate  to  facilitate  its  numerical  solution. 

B.  Statement  of  Problem 


Let  us  consider  the  detonation  transfer  between  two  adjacent  munition 
rounds  as  shown  in  Figure  1  (corresponding  to  an  experimental  configuration 
of  Reference  1).  In  a  gap-test  analogy,  we  may  refer  to  one  round  as  the  donor 
and  the  other  as  the  acceptor.  Of  course,  the  explosives  in  both  the  donor 
and  acceptor  are  essential  for  the  detonation  transfer  here.  Casings,  ambient 
air,  and  any  plastic  shield  form  the  "gap"  or  shock  attenuator.  For  com¬ 
putational  purposes,  we  further  idealize  the  interround  layout  to  two  dimensions 
as  shown  in  Figure  2.  The  symbols  R  ,  R,  r,  and  h  are  defined  in  Figure  2, 
and  we  use  the  parameters  R/R0j  r/R^?  and  h/RQ  to  characterize  our  computational 

problem.  We  shall  only  consider  interround  separations  within  the  range 
1<R/Rq^2,  since  this  is  the  region  in  which  shock  initiation  is  expected  to 

dominate.  Note  that  the  lower  limit  (R/Rq  =  1)  represents  the  mode  of  propa¬ 
gation  of  the  gap  test.  Above  the  upper  limit,  fragment  impact  dominates  the 
initiation  process.  Between  these  two  extremes,  shocks  caused  by  casing  impact 
contribute  to  the  processes  of  interround  communication,  and  we  have  used 
R/Rq  =  1.5  as  a  typical  value  for  sample  calculations.  In  another  group  of 

computations  with  R/Rq  -  1-5,  plastic  shields  of  various  thicknesses,  h/R^, 

are  included  in  the  middle  of  the  air  gap.  We  have  considered  three  basic 
types  of  interround  problem  as  summarized  in  Tables  I  through  III.  These 
consist  of  twenty  cases  all  with  R  =  52.5  mm  and  composition-B  as  the  donor 
charge  ihitiated  on  its  axis.  u 

The  five  cases  of  Table  I  refer  to  two  rounds  in  contact  (R/Rq  =  1*0, 
h/RQ  =  0).  Columns  4-9  describe  the  computational  specifications  for  the 


TABLE  I 

ROUNDS  IN  CONTACT 

AX, AY  IJ  At  N 


Case 

r/R 

o 

Acceptor 

(mm) 

I 

J 

(cells) 

(psec) 

(cycles) 

A-l 

0.15 

PBX-9205 

1.944 

50 

135 

6,750 

0.05 

600 

A- 2 

0.20 

PBX-9205 

1.944 

5Q 

135 

6,750 

0.Q5 

600 

A- 3 

0.27 

PBX-9205 

1.944 

50 

135 

6,750 

0.05 

600 

A- In 

0.15 

Inert 

PBX-9205 

1.944 

50 

135 

6,750 

Q.  05 

700 

A- 3  r. 

0.27 

Inert 

PBX-9205 

1.981 

50 

130 

6,500 

0.05 

900 
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2DE  code,  which  will  be  explained  later.  The  nonreactive  computations  A-ln 
and  A-3n  were  performed  to  examine  shock  loading,  with  all  input  data  the  same 
as  in  A-l  and  A-3  respectively. 

The  ten  computations  of  Table  II  describe  the  second  type  of  sensitivity 
problem,  i.e.,  two  rounds  separated  by  an  air  gap  (R/Rq  =  1.5  and  2.0,  h/RQ  =  0) . 

TABLE  II 

ROUNDS  SEPARATED  BY  AIR  GAP 


Case 

r/Ro 

R/Ro 

Acceptor 

AX, AY 
(ram) 

I 

J 

IJ 

(cell s) 

At 

(ysec) 

N 

(cycles) 

B-l 

0.15 

1.5 

PBX-9205 

1.981 

SO 

150 

7,500 

0.05 

900 

B-2 

0.20 

1.5 

PBX-9205 

1.944 

SO 

150 

7,500 

0.05 

900 

B-3 

0.27 

1.5 

PBX-9205 

1.981 

50 

150 

7,500 

0.05 

900 

C-l 

0.15 

2.0 

PBX-9205 

1.981 

50 

160 

8,000 

0.05 

1,000 

C-2 

0.20 

2.0 

PBX-9205 

1.944 

SO 

165 

8,250 

0.05 

1,000 

C-3 

0.27 

2.0 

PBX-9205 

1.981 

50 

160 

8,000 

0.05 

1,000 

B-ln 

0.15 

1.5 

Inert 

1.981 

SO 

145 

7,250 

0.05 

900 

B-3n 

0.27 

1.5 

Inert 

1.981 

SO 

145 

7,250 

0.05 

900 

C-  In 

0.15 

2.0 

Inert 

1.981 

50 

155 

7,750 

0.05 

900 

C-3n 

0.27 

2.0 

Inert 

1.931 

50 

155 

7,750 

0.05 

900 

Five  cases  with  plastic 

shields 

and  R/  R  = 
o 

1.5  are 

listed 

in  Table  III. 

These 

are  similar 

to  tne  8's 

of  Table  II 

except 

that  a 

plastic 

shield  is  in 

serted  at  the  midpoint  of  the  air  gap. 


TABLE  III 

ROUNDS  SEPARATED  BY  PLASTIC  SHIELD 


Case 

r/Ro 

h/Ro 

Acceptor 

AX,  AY 
(ram) 

I 

J 

IJ 

(cells) 

At 

(ysec) 

N 

(cycles) 

D-l 

0.20 

0.20 

PBX-9205 

1.810 

50 

150 

7,500 

0.04 

800 

D-2 

0.20 

0.28 

PBX-9205 

1.810 

50 

150 

7,500 

0.04 

800 

D-3 

0.20 

0.38 

PBX-9205 

1.810 

50 

150 

7,500 

0.04 

800 

D-  4 

0.20 

0.47 

PBX-9205 

1.810 

50 

150 

7,500 

0.04 

900 

D-5 

0.1S 

0. 14 

PBX-9205 

1.810 

50 

150 

7,500 

0.04 

1,000 
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In  all  the  computations  summarized  in  Tables  I  through  III,  only  five 
materials  are  involved:  composition-B  (64%  RDX,  36%  TNT),  PBX-9205  (92%  RDX, 

6%  polystyrene,  2%  DOP) ,  steel,  air,  and  Plexiglas.  In  order  to  perform  the 
computations,  we  used  the  equation-of-state  and  shock-initiation  constants  of 
Reference  9.  Such  a  large  number  of  material  constants  may  appear  rather  in¬ 
convenient  to  someone  who  is  not  familiar  with  the  computer  code  2DE.  It 
should  be  noted  that  the  constants  have  been  all  calibrated  with  established 
data  from  reliable  sources  (e.g.,  handbooks,  journal  papers,  laboratory  reports), 
In  fact,  several  subprograms  of  2DE  are  executed  one  after  the  other,  each 
calling  for  one  or  two  dozen  of  such  constants  as  input.  Combined  numerical 
results  turn  out  to  offer  a  close  simulation  of  many  detonation  phenomena. 

From  a  computational  point  of  view,  it  is  not  excessive  to  use  the  set  of  con¬ 
stants  as  mentioned. 

II.  NUMERICAL  APPROACH 
A.  Reactive  Hydrodynamic  Code  2DE 

The  computer  code  adopted  for  this  investigation  is  2DE,  a  two-dimensional 
Eulerian  hydrodynamic  code  developed  at  Los  Alamos  National  Laboratory  for 
reactive  flow  problems7’9.  It  is  a  large  code  with  several  features  not  avail¬ 
able  in  other  codes  (e.g.,  HEMP,  HELP).  The  following  are  the  nine  equations 
which  are  used  in  2DE  for  describing  the  reactive  flow  of  detonation  phenomena: 

lP+u3o.*Vi£-=  -  pf~  *  --  ♦  ill  u|  O) 

3t  U  3R  32  PV3R  3Z  R  7 

■'&  (S2Z.  p)  .^SRZ  (3, 

^*^*v>J)  =  .r(|M7Vu)  .sRR(»“  - 

T-1")  (« 

P  =  p  +  q  =  p  (p.  I,  W)  ♦  q  (iCp  Kp  |>-)  (5) 

1  /aw  3 IV  3W  \  l  2  n  \ 

K  \3t  +  U  3R  +  V  3Z  /  exi>  iCo  +  C1P  +  C2P  +-,+  CnP  )  C 6 ) 

3srr  _  /3U  1  ap  \  ^  <.  /  3U  3V  N 

3t  -  2  M  \  VR  +  Tp  3t  )  *  SRZ  {  3Z  -  3R  )  (7) 
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where  p,  U,  V,  W,  P,  I,  S  ,  S_„  and  S  _  denote  the  density,  radial  velocity, 

axial  velocity,  mass  fraction  cf  unreacted  explosive,  pressure,  internal  energy, 
radial  deviatoric  stress,  axial  deviatoric  stress,  and  shear  stress  respectively 
(R,  Z,  and  t  being  the  pertinent  coordinates  and  time  with  a  =  1  for  rectangular 
and  2  for  cylindrical  coordinates).  The  total  pressure  P  is  the  sum  of  two 
parts:  the  thermodynamic  pressure  p  and  the  artificial  viscous  stress  q  for 
the  treatment  of  shock  waves  (with  coefficient  K  =  constant).  Five  subroutines 
are  available7’8  for  the  equation-of-state  computation.  These  are  called  HOM, 
H0M2S,  HOMSG,  H0M2G,  and  H0M2SG.  They  are  used  for  computing  simple  and  mixed 
cells  in  pressure  and/or  temperature  equilibrium  (temperature  T  is  calculated 
according  to  the  Walsh-Christian  scheme).  The  reaction  rate  associated  with 
heterogeneous  shock  initiation  of  an  explosive  is  given  by  Equation  (6)  with 
n  =  14  or  15  (i.e..  Forest  Fire  burn7  using  fourteen  or  fifteen  empirical 
constants).  Note  that  this  burn-rate  equation  is  computationally  coupled 
with  the  equation  of  state  p  =  p  (p,  I,  to').  Elastic-plastic  strength  of  solid 
material  is  taken  into  account  with  von  Mises  yield  criterion  to  compute  SnD, 

and  (note  u  =  shear  modulus).  Other  options  are  also  available, 
including  heat  conduction,  other  bum  models,  and  real  viscosity7. 

Only  two  coordinate  systems  ((X,  Y),  (R,Z))  are  applicable  in  2DE,  namely, 
the  slab  (a  =  1)  and  cylindrical  (a  =  2)  geometry.  Co-planar  computation  (as 
in  our  interround  problem)  can  be  treated  with  circular  input8.  To  set  up  a 
given  problem  for  2DE  computation,  we  first  divide  the  flow  field  into  a 
number  of  rectangular  regions  according  to  the  problem  configuration.  Each 
rectangle  is  again  divided  into  cells,  and  the  total  number  of  cells  (as  IJ  in 
Tables  I- 1 1 1 )  is  determined  by  the  chosen  cell  sizes  (AX  =  X/I,  AY  =  Y/J)  .  The 
time  step  is  At ^ AX/40.  Input  data  (numerical  or  logical)  are  supplied  for 
each  rectangle,  and  the  total  number  of  computing  cycles  is  specified.  Com¬ 
putation  is  then  carried  out  in  six  phases  during  each  cycle  to  yield  a  field 
solution  of  the  flow  problem.  The  finite-difference  scheme  for  all  six  phases 
is  described  in  detail  by  Mader7. 

Let  us  consider  the  numerical  solution  in  the  form  <*>  =  (R,  Z,  t)  with 

$  =  P,  p,  T,  I,  U,  V,  and  to'  (i  =  1,  2,f ••)•  Using  computer  graphics,  we  can 
now  generate  many  plots  to  describe  the  reactive  flow  in  one  (ID),  two  (2D), 
and  three  (3D)  dimensions.  Thus,  we  can  have  3D  plots,  P  =  F^  (R,  Z)  with 

t  =  constant  and  P  =  (Z,  t)  with  R  =  constant.  Isobars  are  2D  plots, 

Fj  (R,  Z)  =  0  with  P,  t  =  constant.  Useful  ID  plots  may  be  P  =  F^  (Z)  with 

R,  t  =  constant,  P  =  F,.  (t)  with  R,  Z  =  constant,  and  P  =  F&  (R)  with  Z,  t  - 

constant.  All  these  plots  are  merely  snapshots  or  histograms.  Likewise, 
computerized  movies  (3D,  2D,  ID)  can  be  produced  by  letting  t  vary  with  R 
and/or  Z. 


B. 


Application  to  NOL  Gap  Test 


We  ran  a  sample  calculation  of  the  NOL  gap  test  prior  to  the  application 
of  the  2DE  code  to  our  numerical  experiments  on  shock  sensitivity  of  munitions. 

A  sketch  of  the  NOL  large  scale  gap  test  (LSGT)  configuration  as  used  in  the  com¬ 
putation  is  shown  in  Figure  3.  In  this  example  VTQ  was  the  test  propellant, 
and  we  set  up  the  problem  with  40  X  140  =  5.6C0  cells.  We  chose  AR  =  AZ  = 

1.826  mm  and  At  =  0.05  usee.  To  account  for  the  material  strength  of  Plexiglas, 
we  also  included  an  elastic-plastic  grid  of  30  X  35  =  1,050  cells.  All 
material  constants  were  taken  from  Reference  9.  The  problem  was  run  in  cylin¬ 
drical  geometry  for  440  cycles,  and  our  computed  results  for  the  shock  initiation 
of  VTQ  with  36.5  mm  gap  are  shown  as  contour  plots  in  Figures  4  and  5.  Note  that 
these  check  with  Bowman's  results9  for  40.2  mm  gap  consistently.  For  com¬ 
plementary  details,  see  Reference  9. 

C .  Application  to  Munitions  Problem 

Since  the  major  events  of  interest  to  us  occur  in  an  almost  coplanar  con¬ 
figuration  as  in  Figure  2,  a  two-dimensional  numerical  simulation  is  justified. 
Thus,  we  ran  our  munitions  problem  in  slab  geometry  with  circular  input  and 
with  elastic-plastic  strength  effects  considered.  In  Tables  I  -  III,  we 
specify  in  the  first  three  or  four  columns  all  types  of  munition  problem  for 
solution,  and  the  remaining  columns  describe  briefly  the  computational  scheme 
with  parameters  already  explained  in  Section  II. A.  It  should  be  remarked 
that  we  modelled  the  detonation  of  the  donor  explosive  (compos it ion-B)  using 
a  C-J  volume  burn  and  that  of  the  acceptor  charge  (PBX-9205)  using  the  Forest 
Fire  burn.  These  models  of  explosive  burn7  differ  in  that  the  former  represents 
a  steady-state  detonation  with  an  infinite  rate  of  energy  release  whereas  the 
latter  is  a  non-steady  shock  buildup  with  a  finite  release  rate,  which  may  or 
may  not  develop  into  a  complete  detonation.  For  the  purpose  of  shock-loading 
analysis,  we  considered  an  inert  acceptor  with  most  of  its  input  data  identical 
to  those  of  a  reactive  acceptor,  except  that  no  Forest  Fire  bum  was  included. 

By  comparing  the  inert  and  reactive  results,  we  can  gain  a  better  understanding 
of  the  mechanism  or  factors  which  control  the  shock  initiation  processes  in 
our  problem.  As  listed  in  Tables  I  -  III,  most  of  the  computations  were  per¬ 
formed  to  treat  detonation  phenomena  with  R/R0  s,  2  and  NAt  50  usee.  These 

values  are  smaller  than  those  associated  with  the  fragmentation  of  the  donor 
cylinder.  Since  the  2DE  code  has  no  provision  to  compute  fracture  of  solid 
material,  we  terminated  our  computations  before  fragmentation  would  have 
occurred. 

Here  a  remark  should  be  made  regarding  the  limitations  of  our  computations. 
The  use  of  slab  geometry  with  circular  input  is  a  practical  approximation  to 
the  real  axi-symmetric  cylindrical  problem.  The  choice  of  casing  thickness 
was  rather  limited,  since  at  least  four  cells  are  needed  to  smear  shock  waves 
with  artificial  viscosity.  A  practical  lower  limit  on  cell  size  is  provided 
by  the  consequent  short  time  step  or  large  number  of  computing  cycles  as  well 
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as  by  the  large  computer  memory  required.  For  these  reasons  we  used  only  r/RQ  = 

0.1S  for  thin-walled  cylinders.  Mixed-cell  computation  of  2DE  code  and  boundary 
irregularities  are  still  far  from  adequate. 

III.  RESULTS  AND  DISCUSSION 
A.  Flow  and  Shock  Patterns  of  Interactive  Munitions 

1.  Reactive  Flows  within  Donor  and  Acceptor. 

Our  twenty  computations  indicate  two  distinct  types  of  reactive  flow 
within  the  donor  and  acceptor.  Once  detonation  is  initiated  from  the  donor  axis, 
diverging  and  converging  cylindrical  shocks  reverberate  in  both  the  gaseous 
detonation  products  and  steel  casing.  The  shock  reverberation  varies  con¬ 
siderably  with  the  parameters  R/RQ,  r/RQ,  an(*  h/R  •  Another  kind  of  Teactive 

flow  is  generated  in  the  acceptor.  This  flow  is  non-cylindrical,  diverging, 
and  shock- induced.  The  threshold  of  shock  initiation  appears  to  be  of  the 
order  of  the  large  scale  gap  test  50%  point.  It  is  the  shock  buildup  that 
shapes  the  contour  of  this  flow.  Shock  loading  is  not  only  complicated  but 
also  different  in  the  two  casings.  Concentric  cylindrical  shocks  reverberate 
in  the  donor  casing,  whereas  only  diverging,  curved  shocks  advance  in  the  acceptor 
casing  from  one  side  to  the  other  -  sometimes  resulting  in  a  collision  which 
sends  backward  running  shock  in  the  explosive.  For  R/RQ  =  1.5  and  2.0,  air 

shocks  are  felt  by  the  acceptor  casing  and  explosive.  However,  they  are  in¬ 
sufficient  to  cause  reaction. 

2.  Propagation  of  Reactive  and  Inert  Shock  Wa  :s. 

The  above  observations  are  common  to  all  the  types  of  problems  wt  have 
studied  (see  Tables  I  -  III).  In  particular,  let  us  refer  to  cases  A-l,  A-ln 
B-l,  B-ln,  C-l ,  and  C-ln.  Figures  6-11  give  their  shock  profiles  plotted 
as  P  =  F2  (X,  t)  along  the  axis  Y  =  0.  These  plots  are  complementary  in  pairs 

(reactive  versus  inert  acceptors).  Thus,  Figures  (>  and  7  show  the  same  diverging- 
converging  flow  on  the  doner  side,  but  the  propagative  responses  differ  on  the 
acceptor  side  depending  on  shock  sensitivity.  In  the  donor  explosive,  a  C-J 
detonation  begins  as  a  spike  of  28  GPa,  spreading  with  a  steady  wave  front  of 
17  GPa  (note  that  cylindrical  divergence  causes  this  reduction  of  11  GPa  in 
the  first  microsecond).  This  steady  detonation  front  advances  with  a  straight- 
line  trajectory  in  the  Xt-plane  (i.e.,  dX/dt  =  detonation  velocity).  About 
5  usee  later  the  detonation  front  is  reflected  from  the  steel  casing,  and  shock 
transmission  and  reverberation  occur  in  the  casing  with  pulses  as  labelled  in 
Figure  7.  A  reflected  shock  in  the  detonation  products  converges  inward  and 
then  diverges  outward,  upon  reflection,  like  a  second  explosion.  (The  converging 
shock  with  a  peak  pressure  of  27  GPa  diminishes  rapidly  at  first  and  then  slowly 
to  1.5  GPa.)  Shock  transmission  to  the  acceptor  (reactive  or  inert)  is  as 
follows.  In  Figure  6,  transmitted  shocks  accelerate  as  soon  as  their  amplitudes 
exceed  the  threshold  of  the  LSGT  50%  point.  The  spread  of  initiation  and 
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Figure  10  -  Shock  Pressure  versus  Distance  and  Time  for  Two  Rounds  Separated  (Reactive 


reaction  can  be  observed  in  Figures  14  and  15  from  t  =  12  usee  or  cycle  240  on. 
The  reactive  shock  in  the  acceptor  builds  up  to  the  C-J  state  (28  GPa)  at  t  = 

16  nsec  when  the  converging  peak  occurs  in  the  donor  gases  (see  Figures  12  and 
14).  We  note  that  this  shock  front  also  describes  a  straight  trajectory  in 
the  Xt-plane  of  Figure  6.  On  the  other  hand,  Figure  7  clearly  shows  the 
penetration  of  shocks  in  the  inert  acceptor  -  stronger  shocks  catching  up 
the  early  transmitted  shocks. 

As  in  Figures  6  and  1 ,  the  reactive  flow  and  shock  waves  on  the  donor  side 
are  the  same  in  Figures  8-11.  In  these  four  problems,  the  air  gap  delays 
the  propagation  of  reaction  or  shock  waves.  Such  prolonged  communication 
furnishes  us  with  a  clear  picture  of:  (a)  air  shocks  picked  up  by  the  acceptor 
casing  and  explosive  as  shown  in  Figure  8,  (b)  different  shock  transmission  to 
the  inert  acceptor  of  Figure  9,  and  (c)  shock  pulses  modulating  in  the  donor 
casing  of  Figures  10  and  11.  Initiation  of  the  acceptor  explosive  in  these 
problems  occurs  as  the  shock  caused  by  casing  impact  builds  to  detonation  as 
shown  in  Figures  8  and  10. 

Computations  A-ln,  B-ln,  and  C-ln  were  performed  in  the  same  way  as  A-l, 
B-l,  and  C-l  except  that  no  reaction  was  permitted  in  the  acceptor.  A  com¬ 
parison  of  these  results  enables  us  to  distinguish  between  the  inert  and  re¬ 
active  responses  of  an  acceptor  to  the  same  shock  stimuli.  Such  a  comparison 
is  useful  to  assess  the  shock  s-nsitivity  of  the  simulated  munition  and  the 
fa.'tors  which  control  the  shock  initiation  of  the  acceptor. 

A  close  look  at  Figures  6-11  indicates  that  early  transmitted  shocks 
cause  no  significant  reaction  in  the  acceptor.  In  all  cases  of  reactive  c>  m- 
munication,  initiation  occurs  in  the  acceptor  as  a  result  of  shock  buildup 
either  by  the  donor's  converging-diverging  flow  with  =  1-0  (see  Figure  12 

for  A-l,  A- 2,  and  A-3  beyond  distance  =  15  era)  or  by  casing  impact  upon  closing 
the  air  gap  with  R/Rq  =  1-5  and  2.0.  Note  that  the  impact  shock  pulses  all 

nave  peaks  exceeding  the  LSGT  50%  point.  As  already  noted,  the  converging- 
diverging  shock  occurs  n  the  donor  as  if  it  were  a  second  explosion.  Figure 
12  shows  the  timing  and  amplitude  of  this  shock  from  a  matrix  of  nine  com¬ 
putations.  We  see  that  both  the  amplitude  and  time  of  occurrence  decrease 
as  the  casing  thickness  increases.  These  will  affect  the  initiation  of  the 
acceptor  explosive.  As  depicted  in  Figure  13,  a  converging  shock  can  also  be 
formed  in  an  inert  acceptor.  Note  that  Figures  12  and  13  are  both  ID  plots 
of  P  =  F4  (X)  along  the  axis  Y  =  0  and  at  specified  times. 

3.  Description  by  Contour  Plots. 

Both  isobar  and  mass-fraction  contour  plots  show  the  reactive  flows 
and  shock  waves  with  greater  clarity.  Without  indicating  the  pressure  levels, 
we  note  that  the  sixteen  isobar  plots  of  Figure  14  show  the  C-J  spike  clearly 
(in  cycle  1),  the  diverging  cylindrical  detonation  front,  reflected  shocks, 
the  converging  shock  (in  cycle  320),  and  the  prolonged  expansion  of  the  de¬ 
tonation  products  and  donor  casing.  The  initiation,  buildup,  and  propagation 
of  detonation  ir  the  acceptor  explosive  from  cycle  240  to  480  are  shown  in 
both  Figures  14  and  IS.  Note  that  here  the  reactive  shock  is  diverging 
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gure  15  -  Mass  Fraction  Contours  of  Two  Rounds 


(although  not  cylindrical)  and  advances  with  decreasing  curvature  (increasing 
strength) .  This  may  be  attributed  to  the  heavy  confinement  on  one  side  and 
the  Taylor  rarefaction  on  the  other.  There  is  appreciable  shock  transmission 
from  the  acceptor  back  to  the  donor.  Note  also  that  an  eccentric  converging 
shock  forms  within  the  acceptor  (see  Figure  14,  cycle  600)  as  a  result  of  the 
detonation  propagation  and  reflection  (see  Figure  14,  cycles  480,  520,  560, 
and  600).  The  irregular  rhape  of  the  donor  casing  in  Figure  15  is  a  combined 
result  of  a  circular  boundary  described  in  rectangular  coordinate  and  mixed¬ 
cell  computations.  Of  course,  this  is  not  the  true  deformation  as  sketched 
in  Figure  16  which  shows  a  superposition  of  two  rounded-off  contour  plots. 

B.  Response  of  Reactive  Acceptors 

1.  Shock  Initiation. 

Shock  waves  penetrating  into  the  acceptor  explosive  compress  it  and 
generate  hot  spots  which  may  multiply  and  cause  significant  reaction.  A  com¬ 
plete  buildup  to  detonation  presumably  depends  upon  either  a  critical  incipient 
energy10  or  certain  reaction  rate.  In  2DE  code  computation  this  rate  is  given 
by  the  Forest  Fire  burn7as  an  exponential  function  of  shock  pressure  only. 

On  the  basis  of  a  single  buildup,  the  Forest  Fire  model  is  constructed  for 
heterogeneous  shock  initiation  as  represented  by  the  Pop  plot  generated  ftom 
wedge-test  results  (note  that  one-dimensional  shock  sensitivity  data  are  used 
here  in  the  same  way  that  1-D  shock  Hugoniots  are  used  in  the  HOM  equation- 
of-state  computations,  i.e.,  without  considering  multidimensional  effects). 
Accordingly,  the  computation  of  the  Forest  Fire  bum  rate  is  an  essential  part 
of  2DE  modeling  of  our  interround  detonation  phenomena.  Figure  17  depicts  a 
heterogeneous  shock  initiation  trajectory  Od  in  the  Xt -plane  as  obtained  from 
a  wedge  test.  The  initial  slope  of  this  trajectory  is  equal  to  the  initial 
shock  velocity  and  the  final  slope  is  equal  to  the  CJ-detonation  velocity 
(curve  od  being  tangent  to  straight  line  dD  along  which  proceeds  steady-state 
detonation).  To  simplify  the  description,  let  us  consider  chord  od  with 
Ua  =  tj,  the  average  velocity  of  shock  buildup  (see  the  insert  in  Figure 

17),  where  is  the  distance  to  detonation  and  t^  the  time  to  detonation. 

The  parameter  Ua  increases  with  the  strength  of  the  shock  stimulus  or  the 

sensitivity  of  the  explosive  (U  D,  D  being  CJ-detonation  velocity).  Con- 

versely,  tardy  responses  correspond  to  small  values  of  Ua  as  in  acceptors 

with  thick  casings  or  at  large  distances  of  separation  (see  Figures  24  and  25). 

2.  Effects  of  Parameter  Variations. 

Shock  initiation  of  the  acceptor  explosive  may  be  influenced  by  many 
factors.  Here  we  consider  those  parameters  which  are  obtainable  from  our  2DE 
computations.  In  Table  IV  we  list  values  of  such  parameters,  including  our 
notations  and  symbols. 


10ii lalkev,  F.  E.  and  Waaley,  R.  J.,  "Critical  Energy  for  Shock  Initiation 
Explooiveatoffe  17,  9  (1969).  See  also  the  Sixth  Syrni^oaium  (International) 
cn  Detonation ,  ACR-221  (1976)  for  several  diBcueeiona  on  this  topic. 
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SHOCK  INITIATION  RESULTS 
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With  D  and  P^j  denoting  the  CJ-detonation  velocity  and  pressure 

respectively,  we  also  introduce  tQ  =  Rq/D  as  a  characteristic  time  for  the 

reactive  acceptor.  Now  we  can  normalize  all  initiation  parameters  of 
Table  IV  and  plot  the  results  as  shown  in  Figure  18  -  25.  Each  of  these 
plots  may  be  looked  upon  as  an  expression  of  the  form: 


ni  =  fi  (R/Ro*  r/Ro‘  h/Ro’  *  5  (10) 

where  (i  =  1,  2,,.,)  are  response  parameters,  say,  =  t^/t^,  =  tcj/t0>* 

and  £;  is  an  explosive  parameter  implicitly  held  constant  here.  We  are  in¬ 
terested  in  the  general  trends  exhibited  in  these  plots  rather  than  in  their 
analytic  expression  such  as  Equation  (10). 

It  is  worth  noting  that  all  curves  in  Figures  18  and  19  have  a  trend 

of  increasing  reaction  delay  as  the  casing  thickness  or  interround  separation 

is  increased.  Specifically,  we  observe  in  Figure  19  greater  reaction  delay 

with  r/.  >  0.2  than  with  r/R  <  0.2  for  R/R  =  1.0.  In  Figure  18  we  note 

o  oo 

little  delay  with  1.0  <  R/R0  <  1.5  but  considerably  greater  delay  with 

1.5  <  R/R  <  2.0  for  r/R  =  0.1S. 
o  o 

in  Figures  20  and  21  a  similar  trend  is  shown  using  the  time  required 
for  shock  buildup  to  detonation.  A  comparison  between  Figures  18  and  20 
reveals  that  for  r/R^  =  0,2  a  change  of  R/Rq  from  1.0  to  1.5  results  in  a 

reaction  delay  (Figure  18)  but  an  advance  of  full  detonation  (Figure  20).  It 
is  also  interesting  to  note  that  in  Figure  21  the  difference  between  R/RQ  = 

1.0  and  1.5  is  narrow  down  with  0.15  <  r/R  <  0.20.  A  similar  trend  exists 

o 

between  R/R  =  1.0  and  2.0  with  0.20  <  r/R  <  0.25. 
o  o 

The  initiation  parameter  X^  was  taken  as  the  distance  traversed  by 

the  reactive  shock  front  along  the  symmetry  axis  during  buildup  to  detonation. 
It  was  then  used  as  a  simple  one-dimensional  model  for  computing  Ua/D  = 

X ,/Dt .  and  P./P_T  =  (X  eA)1//fi/p_T  for  Figures  22  -  25,  where  A  and  B  are 
da  .1  UJ  d  /  LJ 

the  Pop  plot  constants  as  taken  from  Reference  9.  Of  course,  the  shock 
stimuli  and  buildup  in  our  interround  problems  are  multiple,  two-dimensional 
waves.  Here  we  seek  only  a  crude  representation  of  the  initiation  phenomena 
using  Figures  22  -  25  (more  realistic  model  not  yet  available).  Thus,  the 
trend  of  Figure  22  may  be  interpreted  as  follows.  The  shock  stimuli  are 
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Figure  21  -  Effects  of  Casing  Thickness  on  Time  to  Detonation 
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Figure  24  -  Effects  of  Interround  Separation  on 
Average  Velocity  of  Shock  Buildup 


due  to  the  converging-diverging  cylindrical  flow  as  if  it  were  a  second  explo¬ 
sion  in  the  donor  with  R/Rq  =  1-0,  but  they  derive  from  quite  different  mecli- 

ansiras  for  R/Rq  =  1.5  and  2.0.  Upon  closing  the  air  gap,  casing  impact  deliver 

to  the  acceptor  explosive  shocks  which  coalesce  as  if  they  were  stronger  initial 
pulses  for  R/Rq  *  1.5.  As  R/Rq  increases  from  1.5  to  2.0,  the  shock  stimuli 

become  weaker  again.  It  is  worth  noting  that  the  change  in  trend  of  Figure  22 
matches  that  of  Figure  20  consistently  -  more  delay  of  detonation  due  to  weaker 
shock  stimuli.  Within  the  range  considered,  there  is  a  minimum  delay  due  to 
the  maximum  stimuli  at  R/Rq  *1.5.  In  Figure  23,  a  decreasing  trend  is  observed 

of  shock  intensity  versus  casing  thickness.  Note  the  relative  shock  intensities 
for  R/Ro  =  1.5,  2.0,  and  1.0  at  r/RQ  =  0.20.  In  fact,  these  are  due  to  different 

mechanisms  of  shock  transmission  as  mentioned  earlier.  Both  Figures  24  and  25 
indicate  that  shock  buildup  is  retarded  either  by  interround  separation  or  by 
casing  thickness,  as  expected  from  intuition. 

Using  Figures  18  -  25,  we  have  examined  the  effects  of  varying  interround 
separation  and  casing  thickness  on  shock  sensitivity  of  munitions.  A  simple 
assessment  is  to  determine  the  general  trends  of  these  figures  which  describe 
relations  between  dimensionless  parameters.  Such  trends  can  serve  as  a  practical 
guide  for  scaling  or  interpolating  detonation  phenomena  within  the  range  of  our 
numerical  simulation. 

C.  Response  of  Inert  Acceptors 

1.  Shock  Loading. 

To  check  the  shock  loading  on  an  acceptor  explosive,  we  computed  the 
dynamic  response  of  the  acceptor  filled  with  a  dummy  explosive.  Shock  pres¬ 
sures  were  taken  for  the  explosive  cell  nearest  the  casing  as  a  function  of 
time  only  -  viz.  P  «  F  (t)  as  mentioned  in  Section  II. A.  P*  refers  to  the 

peak  pressure  at  time  t„  of  the  transmitted  shock  pulse;  P  is  the  peak  pressure 

at  time  t  of  the  impact  shock  pulse;  p'  and  p'  are  shock  pressures  at  times 

l  /  2 

tj  and  respectively  as  mentioned  earlier  (see  Table  V  );  and  P  dt  is  integrated 
over  time  in  order  to  apply  the  critical  energy  criterion10.  Note  the  dif¬ 
ference  between  and  P„  with  <  t„.  These  are  usually  too  weak  to  induce 
any  significant  reaction.  With  t  < 't',  p'  belongs  to  the  impact  shock  pulse 

^  ^  /X 

which  is  responsible  for  shock  initiation  whenever  P  exceeds  the  LSGT  50% 

point.  Notably,  as  calculated  for  the  inert  acceptor,  is  more  comparable 

with  P^  than  with  the  LSGT  50%  point.  But  P2  values  for  the  reactive  acceptor 

turn  out  to  be  just  about  at  the  LSGT  50%  point  (see  Table  IV).  Such  com¬ 
parison  is  consistent  with  the  shock  buildup.  Table  V  gives  the  numerical 
results  of  inert-acceptor  response,  with  all  t's  in  ysec  and  all  P's  in  GPa. 
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TABLE  V. 


SHOCK  LOADING  RESULTS 


Case 

F/R0 

r/Ro 

t* 

P* 

/\ 

t 

P 

m* 

/P2dt 

*1 

1 

P1 

t2 

» 

P2 

A- In 

1.0 

0.15 

10.2 

4.0 

19.4 

7.9 

6.1 

9 

1.3 

12 

2.3 

B-ln 

1.5 

0.15 

23.0 

0.8 

35.8 

6.9 

7.5 

17 

0.07 

32 

1.2 

C-ln 

2.0 

0.15 

25.2 

0.1 

43.6 

5.5 

12.0 

25 

0.1 

42 

1.0 

A-3n 

1.0 

0.27 

13,2 

1.1 

24.2 

8.6 

15.0 

11 

0.1 

20 

1.3 

B-3n 

1.5 

0.27 

20.2 

0.1 

36.2 

7.9 

4.0 

19 

0.04 

- 

- 

C-3n 

2.0 

0.27 

23.8 

0.07 

38.6 

2.3 

2.5 

28 

0.05 

- 

- 

It  is  apparent  that  both  the  inert  and  reactive  acceptors  should  re¬ 
spond  identically  to  equal  shock  stimuli  which  cause  no  initiation  at  all 
(see  Figures  10  and  11  and  lines  C-l  and  C-ln  in  Tables  IV  and  V  respectively). 
Any  differences  between  values  of  Tables  IV  and  V  are  due  to  chemical  reaction 
which  affects  shock  buildup,  and  we  can  check  them  with  shock  sensitivity  data. 
In  fact,  the  response  of  an  inert  acceptor  involves  a  complicated  pattern  of 
shock  propagation  and  interaction.  Here  we  have  simplified  our  analysis  of 

shock  loading  by  considering  only  P„  and  P  which  are  essential  for  detonation 
propagation.  Below  we  seek  to  examine  our  shock  loading  results  from  a  view¬ 
point  which  is  complementary  to  Section  III.B.2. 

2.  Effects  of  Parameter  Variations. 

Using  the  numerical  results  presented  in  Table  V,  we  can  now  examine 
the  shock  behavior  of  inert  acceptors  in  order  to  gain  some  insight  into  the 
problem  of  detonation  transfer.  To  construct  the  plots  of  Figures  26  -  28, 
we  have  norsialized  the  pertinent  parameters  with  the  following  standard 
constants5’11  for  PBX-920S: 

PCj  =  28.1  GPa,  D  =  8.13  ram/usec  (Y  =  2.97) 

X  *  50.8  mm,  Px  *  5.0  GPa  (LSGT  50%  point) 

P2t  =  5,100  (GPa)?psec 

LJ  o  2 

E  a  0.0019  (normalized  P  t  value) 


11 Dobratz ,  B.  M,  (ed). 
Simulants , "  Laurence 


"Properties  of  Chemical  Explosives  and  Explosive 
Livermore  Laboratory  Report  UCRL-S2997  (March  1981). 
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2  2 

The  last  constant  is  obtained  from  E  =  P  t/PrTt  =  Uu/Du  with  U  =  3.9  mm/ysec 

wj  o  2 

(using  P^  and  shock  Hugoniot  to  compute),  w  =  Put  =  P  t/pQU  =  14.6  GPa-nun 

(using  the  value11  for  composition-B  as  a  close  estimate  here),  and  = 

P.TR  /(Y+l)  =  372  GPa-mm. 
uJ  o 

Figure  26  shows  that  all  shock  pressures  P#  and  ‘p'  decrease  as  the 
interround  separation  R/Rq  increases.  Using  the  LSGT  50%  point  (viz.  horizontal 

line  P  /P  =  0.18)  as  a  guide,  we  see  that  the  transmitted  pulses  P*  are  too 

X  CJ 

weak  and  the  impact  shock  pulses  P  are  responsible  for  detonation  to  occur  in 
reactive  acceptors.  These  can  readily  be  verified  by  comparing  Figures  6  and 
8  with  Figures  7  and  9  respectively.  Without  including  P*/P„T  which  is  not 

/\ 

critical  here.  Figure  27  describes  the  variation  of  P/P„T  with  r/R  .  We 

LJ  O 

note  first  an  increasing  trend  for  R/R  =  2.0.  With  r/R  increased  from 

o  o 

0.1S  to  0.27,  the  shock  sensitivity  of  a  reactive  acceptor  at  R/RQ  =  2.0  will 
be  reduced  to  a  "no-go"  level  (see  Figures  26  and  27). 

The  critical  energy  criterion10  has  often  been  used  as  a  guide  for 
the  assessment  of  shock  initiation.  The  particular  constant  is  the  shock 

2 

energy  per  unit  area,  u>  =  Put  =  P  t/pQU,  due  to  plate  impact  in  sensitivity 

test.  Thus,  we  arrived  at  the  value  E  =  0.0019  earlier.  Using  this  and  values 
of  Table  V,  we  plot  the  horizontal  (dotted)  line  and  curves  in  Figure  28  as  a 

check  of  shock  sensitivity  according  to  the  P2t  criterion.  Thus,  the  horizontal 
line  may  be  regarded  as  the  initiation  threshold  f  rove  which  we  would  predict 
"go"  and  below  which  "no  go."  A  comparison  bet’  \  such  prediction  and  our 

results  of  Table  IV  indicates  three  points  of  Fi&-re  28  (viz.  a'  ,  c'  ,  c)  in 

f 

agreement  and  others  (viz.  a,  b,  b  )  not  so.  Noti  that  the  separation  R/R  = 

2  ° 

1.5  will  provide  the  smallest  chance  of  initiation,  if  a  P  t  criterion  is 
applicable.  However,  as  noted  in  Table  IV,  initiation  occurred  for  R/Rq  = 

1.5  as  well,  and  the  P2t  criterion  may  not  apply  here. 

D.  Effects  of  Plastic  Shield  between  Rounds 

It  has  been  shown  elsewhere1'5  that  plastic  shields  between  rounds  were 
able  to  suppress  detonation  transfer,  provided  their  thickness  and  spacing 
were  adequate.  We  pursued  a  numer"  %al  study  with  values  of  the  shield  thickness, 
h  =  10,  IS,  20  and  25  mm,  for  R/R  =  1.5  and  r/RQ  =  0.20  (see  Table  III).  Our 

numerical  results  are  given  in  Table  VI,  including  D-0  which  is  in  fact  B-2  of 
Table  IV,  and  D-5  for  the  thinnest  shield  considered  here.  We  have  confined 
our  attention  to  R/R  =  1.5  and  r/R  =  0.20  (see  Figures  22  and  24  for  the 
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Figure  26  -  Effects  of  Interround  Separation  on  Inert  Shock  Loading 


Ill 


reduction  of  P.  and  U  respectively).  In  Figure  29  is  plotted  the  only  complete 

X  Si 

set  of  data  on  a  normalized  basis.  Unfortunately,  almost  all  computations 
terminated  early.  With  six  data  points,  we  obtained  an  N-shaped  curve  describing 
vhe  effect  of  increasing  shield  thickness  from  one  extreme  value  to  the  other. 

It  is  interesting  to  note  the  increase  of  t  /tQ  Figures  29  and  18.  With 

r/RQ  =  0.20  (and  h/R^  =  0),  a  40%  increase  is  realized  in  Figure  18  by  changing  R/RQ 

from  1.5  to  2.0.  With  r/R  =  0.20  and  R/R  =  1.5,  a  160%  increase  is  shown  in 

o  o 

Figure  29  by  changing  h/RQ  from  0  to  0.20.  Although  this  part  of  our  study 

is  far  from  conclusive,  the  available  results  indicate  remarkable  benefit  of 
thin  shields. 

IV.  SUMMARY  AND  CONCLUSIONS 

This  report  presents  the  results  of  our  exploratory  numerical  experiments 
on  shock  sensitivity  of  munitions.  We  have  used  the  reactive  hydrodynamic 
code  2DE  to  evaluate  the  posed  problem12  which  covers  three  modes  of  detonation 
transfer,  with  twenty  computations  all  listed  in  Tables  I  -  III.  The  qualif¬ 
ication  of  2DE  as  a  tool  for  numerical  simulation  was  discussed  at  the  outset. 

Our  twenty  computations  describe  the  munition  problem  in  de  ail  which  is  not 
readily  accessible  by  other  approaches.  Two  distinct  reactive  flows  are  in¬ 
volved  in  the  analogue,  with  shock  waves  communicable  between  them.  The 
numerical  aspects  of  these  can  serve  to  explain  the  mechanism  and  development 
of  the  detonation  phenomena. 

Shock  waves  play  an  essential  role  in  all  phases  of  the  flow  and  reaction. 

Inert  shocks  are  propagated  in  the  gaseous  products,  steel  casings,  ambient 
air,  and  plastic  shield  -  with  patterns  complicated  by  various  interfaces  or 
boundaries.  Thus,  cylindrical  shocks  behave  distinguishably  within  the  donor; 
their  converging  and  diverging  in  the  gases  resemble  a  second  explosion;  and 
their  reverberation  in  the  casing  spreads  fluence  forward  and  backward.  Upon 
closing  the  air  gap  (if  any),  casing  impact  add  multiple  shocks  to  the  scene. 

One  after  the  other,  shock  pulses  arrive  at  the  acceptor  explosive;  it  is 
the  impact  shock  pulse,  with  a  peak  exceeding  the  LSGT  50%  point,  that  evolves 
reactively  into  CJ-detonation,  for  the  configurations  investigated. 

Our  numerical  results  further  offer  a  parameter  analysis  for  shock  initiation 
and  loading,  with 

■  *i  <R/V  r/V  h/VE  J 

depicted  in  Figures  18-29.  A  general  trend  of  Figures  18  -  25  indicates 
that  shock-sensitivity  parameters  n  (i  =  1,2-..)  tend  to  decrease  as  R/Rq  or 

r/R  increases.  Such  a  trend  may  be  used  for  practical  purposes  to  scale  or 
o 


lzHowe}  P.  M. ,  Huang ,  Y.  K. ,  and  Arbuokle ,  A.  L.,  "A  Numerical  Study  of  De¬ 
tonation  Propagation  between  Munitions , "  The  Seventh  Symposium  (International ) 
on  Detonation }  Vol  I  (Preprint) 3  pp  512-517  (June  1981). 
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to  interpolate  detonation  phenomena  within  the  range  of  our  numerical  simulation, 
Figures  26  and  27  highlight  the  hypothesis  that  the  impact  shock  pulse  is  essential 
for  shock  initiation.  Yet  it  calls  for  special  attention  to  check  Figure  26  with 
Figure  22.  At  the  first  sight  they  seem  to  contradict  each  other.  However,  P. 

/\  1 
and  P  refer  to  shock  stimuli  in  two  different  phases  of  the  initiation  process, 

and  one  is  bound  to  top  the  other.  It  is  understandable  that  an  out-of-phase 

2 

trend  prevails  during  shock  buildup.  The  P  t  criterion  is  shown  in  Figure  28  to 

indicate  the  least  chance  for  initiation  to  occur  at  R/Rq  =  1.5,  but  it  fails  to 

confirm  our  result  as  validated  with  Figures  18  -  27.  The  trend  of  Figure  29 
shows  the  remarkable  benefits  of  chin  plastic  shields.  It  should  be  remarked 
that  shock  initiation  can  be  affected  by  many  factors  -  chemical,  hydrodynamic, 
and  structural.  Our  analysis  has  placed  more  emphasis  on  structural  parameters, 
and  further  extension  will  have  to  be  made. 

For  lack  of  comparable  data,  the  accuracy  of  our  numerical  results  may  be 
questioned;  but  there  should  be  no  doubt  about  the  state  of  art  for  our  modeling7. 
It  is  anticipated  that  such  situation  will  be  improved  in  our  subsequent  work. 
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